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ABSTRACT 

The  external  chemical  reactivity  of  graphene  sheet,  fullerenes  and  carbon  nanotubes  has  been 
investigated.  The  total  reaction  energy  is  analyzed  with  several  contributing  terms  and 
formulated  as  a function  of  the  pyramidal  angles  of  C atoms.  We  have  determined  the  parameters 
for  the  formulae  from  ab  initio  simulation  of  graphene.  We  have  applied  them  to  predict 
hydrogenation  energy  of  several  nanotubes  and  C6o,  and  demonstrated  that  the  predicted  total 
reaction  energies  are  very  close  to  the  results  of  total  energy  pseudo-potential  density  functional 
theory  calculations.  This  analysis  can  be  used  to  predict  the  reaction  energy  and  local  bonding 
configuration  of  a reactant  with  diverse  fullerenes  and  nanotubes  within  0.1  eV  accuracy. 


INTRODUCTION 

There  has  been  much  research  interest  in  carbon  nanotubes  and  fullerenes  since  the  discovery 
of  C6o  [1].  They  have  been  considered  as  promising  materials  for  nanotechnology  applications, 
such  as  biochemical  and  gas  sensors  [2]  and  molecular  transistor  [3].  From  the  recent  studies  on 
possible  nanodevice  applications,  it  has  been  recognized  that  the  surface  functionalization  of 
nanotubes  and  fullerenes  would  play  an  important  role  for  nanodevice  development.  In  order  to 
functionalize  nanotubes  and  fullerenes,  the  chemical  reactivity  of  carbon  atoms  need  to  be 
understood  with  a quantitative  accuracy. 

Generally,  the  chemical  reactivity  on  the  external  surface  of  a fullerene  or  a nanotube  is 
characterized  by  local  bonding  configuration  of  carbon  atoms,  more  specifically, 
pyramidalization  (0P)  of  C atoms  as  illustrated  in  Fig.  1 [4,  5],  Since  the  surface  of  a fullerene  or 
a nanotube  is  curved,  it  is  natural  to  have  pyramidalized  C atoms  as  shown  in  Fig.  1 . 
Pyramidalization  changes  the  hybridization  of  atomic  orbitals  at  the  C atom  so  that  the  7t  orbital 
contains  different  portion  of  s and  p orbitals  leading  to  different  chemical  reactivity.  For 
example,  graphite  has  planar  structure  (&p  = 0)  corresponding  to  sp2  for  a bonds  and p for  7t 
bond.  However,  fullerenes  and  nanotubes  have  the  hybrid  bonding  orbitals  between  sp 2 to  sp 3. 
Because  of  this  hybridization,  fullerenes  and  nanotubes  are  known  to  be  more  reactive  than 
graphite.  The  degree  of  pyramidalization  ( 0P)  is  defined  by  the  angle  between  o bond  and  n 
orbital  and  named  as  pyramidal  angle  (Fig.  1). 

When  a C atom  interacts  with  an  external  chemical  reactant,  several  processes  occur 
simultaneously:  orbital  hybrid  changing  toward  sp3,  n bond  breaking,  and  reaction  between  free 
n orbital  and  external  reactant.  In  order  to  quantify  the  analysis  of  the  chemical  reactivity,  we 
divide  a reaction  into  several  contributing  parts:  (a)  straining  of  surface  C atom,  which  is 
changing  the  hybrid  of  atomic  orbital  of  C atom,  (b)  breaking  of  Jt  bond  and  binding  with 
external  reactant,  and  (c)  local  relaxation  of  neighbor  C atoms.  We  name  several  energy  terms 
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following  each  contributing  part,  such  as,  "strain  energy"  for  (a),  "binding  energy"  for  (b),  and 
“local  relaxation  energy”  for  (c).  To  elucidate  these  terms,  we  will  use  a hydrogen  atom  as  a 
point  probe  to  measure  binding  energy.  We  have  analyzed  and  formulated  these  energies  in 
terms  of  pyramidal  angles.  And  we  have  applied  them  to  graphene  sheet,  several  nanotubes,  and 
C6o  fiillerene.  We  have  chosen  graphene  sheet  as  a representative  system  for  large  fullerenes  and 
nanotubes.  We  have  chosen  several  nanotubes  with  different  radii  such  that  the  pyramidal  angles 
are  between  those  of  graphene  sheet  and  Cm-  In  the  analysis,  we  have  calculated  the  energy 
parameters  for  our  formulae  from  graphene  sheet  data.  These  parameters  and  the  energy 
expressions  are  used  to  predict  the  interaction  energies  of  nanotubes  and  C6o-  The  accuracy  of  the 
prediction  is  tested  by  comparing  to  full  ab  initio  simulations. 


Figure  1.  Pyramidal  angle  ( 6P ) is  defined  by  the  angle  between  7t  orbital  and  o bond  minus  90°. 


COMPUTATIONAL  DETAILS 

All  the  results  in  this  work  are  calculated  using  total  energy  pseudo-potential  density 
functional  theory  (DFT)  method  [6,  7],  and  a supercell  approximation  is  used  to  simulate  the 
periodic  systems.  In  order  to  minimize  the  interactions  between  the  repeated  images,  we  include 
5 A vacuum  separating  the  neighboring  images.  Kohn-Sham  single-electron  wave  functions  are 
expanded  by  50,000  plane  waves  for  C6o,  and  between  1 0,000  and  36,000  for  nanotubes 
depending  on  their  diameters,  and  1 8,000  for  graphite,  corresponding  to  40  Ry  cut-off  energy. 
We  use  single  &-point  for  C6o,  because  it  is  a molecule.  For  a graphene  sheet  and  nanotubes, 
i-point  convergence  was  tested  and  it  is  confirmed  that  a single  £-point  sampling  for  graphene 
sheet  and  four  ft-point  for  nanotubes  are  enough  for  the  convergence.  During  geometry 
optimizations,  the  atomic  positions  were  relaxed  until  the  forces  became  smaller  than  0.05  eV/A. 

The  initial  values  of  the  pyramidal  angles  do  of  graphene  sheet,  C nanotubes,  and  C6o  used  in 
this  work  are  listed  in  Table  I.  When  we  compute  the  strain  energy,  we  fix  all  the  C atoms  and 
move  one  C atom,  on  which  the  reaction  will  occur,  normal  to  the  surface  in  order  to  increase 
pyramidal  angle  artificially.  We  calculate  the  total  energy  of  the  strained  structure  and  get  the 
strain  energy  of  given  structure  from  the  difference  of  total  energies  between  the  strained 
structure  and  the  initial  structure.  In  order  to  understand  strain  energy  dependency  on  pyramidal 
angle,  we  have  computed  a set  of  strain  energy  with  different  pyramidal  angles. 

We  have  calculated  the  binding  energy  of  H atom  for  each  strained  structure.  We  have 
introduced  an  H atom  above  the  strained  C atom  structure  generated  in  the  previous  analysis  and 
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calculated  hydrogenation  energy  to  compare  the  chemical  reactivity  difference  induced  by 
pyramidal  angle  change.  During  the  calculation  of  binding  energies,  we  keep  the  C atoms  fixed 
to  maintain  the  artificially  produced  pyramidal  angle  and  atomic  hybrid,  and  relax  the  position  of 
the  introduced  H atom.  We  have  also  simulated  fully  relaxed  hydrogenation  energy  to  compare 
with  our  predictions.  We  relax  C atoms  up  to  the  third  nearest  neighbors  for  graphene  sheet,  up 
to  the  first  nearest  neighbors  for  nanotubes,  and  all  the  C atoms  for  C6o  fullerene  during  the  full 
relaxation  simulations. 


Table  I.  Initial  Pyramidal  Angles  (9o)  of  graphene  sheet,  several  nanotubes,  and  C6o  (rad.) 


Oo 

9o 

Graphite 

0.000 

(10,0) 

Nanotube 

0.090 

(15,0) 

Nanotube 

0.059 

(8,0) 

Nanotube 

0.113 

(12,0) 

Nanotube 

0.075 

C60 

0.202 

RESULTS  and  DISCUSSION 


When  we  strain  one  C atom,  the  pyramidal  angle  and  the  bond  length  of  the  pulled  C atom 
are  changed.  So,  we  have  divided  strain  energy  into  two  terms,  bending  energy  and  bond 
stretching  energy.  We  have  analyzed  those  terms  by  force  constant  models  and  described  as 
quadratic  function  of  pyramidal  angle  ( 9P ) changes  from  the  initial  equilibrium  pyramidal  angle 
(Q0)  and  bond  length  (a)  changes  from  the  initial  equilibrium  bond  length  (ao),  respectively. 
When  C atom  is  pyramidalized  as  shown  in  Fig.  1 , the  bottom  triangle  is  consist  of  neighboring 
three  C atoms.  Because  the  triangle  of  bottom  pyramid  is  not  changed  during  the  straining  and 
the  strained  C atom  is  pulled  up  normal  to  the  triangle,  a cos (9P)  remains  a constant  and  bond 
length  (a)  can  be  replaced  by  a function  of  pyramidal  angle  (6P)  as  shown  in  equation  (1).  Due 
to  three  bonds  in  strained  C atom,  there  are  factors  of  three  in  equation  (1).  Since  nanotubes  and 
fullerenes  are  graphite-like  materials,  we  develop  the  scheme  for  graphene  sheet  and  apply  it  for 
the  other  carbon  materials  with  averaged  bond  lengths  and  averaged  pyramidal  angles  due  to  Cj 
symmetry  in  graphene  sheet. 


E strain  ~ ^ ^ bending  (Pp  @0  ^ + 2^ ‘ 


slretching^O 


rcos{90)  A 


-1 
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DFT  simulation  of  strain  energy  for  graphene  sheet,  (10,  0)  nanotube,  and  C6o  is  shown  in 
Fig.  2.  We  fit  the  DFT  simulation  data  for  graphene  sheet  to  equation  (1)  and  determine  the 
spring  constants:  10.70  eV/rad2  for  handing  and  7.44  eV/A2  for  kstretchmS ■ In  Figure  2,  we  plot  the 
predicted  strain  energy  curves  for  (10,0)  nanotube  and  C6o  and  these  curves  show  a good 
agreement  with  the  DFT  simulation  results. 
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Figure  2.  DFT  simulation  results  for  strain  energy  and  binding  energy.  Strain  energy  is  always 
positive  and  shown  as  top  three  set  of  data  and  binding  energy  is  shown  as  bottom  three  set  of 
data.  We  use  square  for  graphene  sheet,  triangle  for  (10,0)  nanotube,  and  circle  for  C6o-  The 
solid  line  for  graphene  is  a fitting  to  the  equations,  and  the  solid  lines  for  the  nanotube  and  C6o 
are  predictions  from  the  equations.  The  dotted  line  is  the  sum  of  strain  energy  and  binding 
energy  for  graphene  sheet. 


Binding  energy  is  composed  of  two  different  contributions,  it  bond  breaking  and  reaction 
between  external  reactant  and  free  7t  state.  We  use  it  bond  breaking  energy  of  graphene  sheet  as  a 
reference  and  represent  it  bond  breaking  energy  of  other  material  based  on  it  bond  breaking 
energy  of  graphene  sheet  and  the  difference,  which  is  shown  in  equation  (2),  where  Em/i  is  the 
energy  difference. 

En  - Egraphtte  + Eshtfl  (2) 

For  the  reaction  of  free  it  state,  we  use  POVA  (n  Orbital  Vector  Analysis)  [8],  which  assigns 
one  s and  three  p orbitals  to  three  G bonds  and  determine  it  state  from  the  rest.  Equation  (3) 
shows  it  state  in  terms  of  the  s and  p orbitals  derived  from  POVA.  When  an  external  reactant  is 
introduced  to  the  7t  state  expressed  in  equation  (3),  the  energy  for  the  reaction  with  an  external 
reactant  (x)  will  be  expressed  by  <s\Hlx>  and  <p\H\x>  using  proper  Hamiltonian  H.  The  binding 
energy  is  the  sum  of  it  bond  breaking  energy  in  equation  (2)  and  the  energy  for  the  reaction  of 
free  7t  state.  And  the  binding  energy  will  be  expressed  as  equation  (4),  where  Esx  and  Epx  are 
<s|//|x>  and  <p\H]x>,  respectively. 


|*)  = & tan(6>»  + A/l  - 2 tan2(0p)|  p)  (3) 

E binding  = v/2  tan (0p  )ESX  + ^1-2  tan  (0p)Epx  + Egraphill  + E 'sh!j,  (4) 
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DFT  simulation  results  of  binding  energy  are  also  summarized  in  Fig.  2.  We  fit  the  graphene 
sheet  data  to  equation  (4)  to  get  the  parameters:  -5.24  eV,  -1 .55  eV,  and  0.87  eV  for  Esx,  Epx  and 
E graphite,  respectively.  The  fitted  binding  energy  curve  is  drawn  in  Fig.  2 as  solid  curve 
connecting  square  data  points.  From  a comparison  between  this  curve  and  other  DFT  data  in  Fig. 
2,  we  can  see  that  a small  energy  shift,  around  -0.2  eV  — 0.3  eV,  is  needed  for  (10,0)  nanotube 
and  C6o  to  match  the  binding  energy  curve  of  graphene  sheet.  In  equation  (4),  the  first  two  terms 
depends  only  on  the  pyramidal  angle  (8P),  and  the  third  term  is  a constant.  The  only  possible 
energy  difference  for  different  structured  materials  is  coming  from  the  last  term,  Eshift-  As  shown 
in  Fig.  2,  Eshf,  is  not  affected  very  much  by  the  initial  pyramidal  angle,  and  we  choose  -0.2  eV  as 
EShift  for  other  nanotube  systems. 

Total  reaction  energy  is  obtained  as  the  sum  of  strain  energy,  binding  energy,  and  local 
relaxation  energy,  as  describe  above.  First,  we  try  to  estimate  the  total  reaction  energy  using  only 
strain  energy  and  binding  energy,  which  are  represented  by  equations  (1)  and  (4)  respectively.  As 
an  example,  we  determine  the  total  reaction  energy  of  graphene  sheet  from  the  minimum  of  the 
sum  of  strain  energy  and  binding  energy  from  Fig.  1.  We  compare  the  total  reaction  energy  from 
full  relaxed  result  and  get  0.16  eV  of  energy  difference,  which  is  the  local  relaxation  energy. 

Finally,  we  now  predict  the  total  reaction  energy  of  an  H atom  on  different  nanotubes  and 
Cots-  We  approximate  that  the  local  relaxation  energy  does  not  strongly  depend  on  the  initial 
structure  and  use  the  value  of  graphene  sheet,  -0.16  eV,  for  other  nanotubes  and  C6o-  Therefore, 
we  determine  the  total  reaction  energy,  especially  hydrogenation  energy  for  this  case,  from 
equation  (1 ) for  strain  energy,  equation  (4)  for  binding  energy,  and  -0. 1 6 eV  for  local  relaxation 
energy.  The  predicted  total  reaction  energy  results  are  summarized  in  Fig.  3 as  square  date 
points.  In  Figure  3,  we  compare  the  predicted  data  with  the  DFT  calculations  with  full  relaxation 
(circle  data  points).  Figure  3 shows  that  we  can  predict  hydrogenation  energy  of  fullerenes  and 
nanotubes  within  0.1  eV  errors  compared  to  the  full  DFT  simulations. 


CONCLUSIONS 

We  have  shown  that  the  total  reaction  energy  can  be  divided  into  three  energy  terms,  strain 
energy,  binding  energy,  and  local  relaxation  energy.  We  also  have  shown  that  those  energy  terms 
can  be  formulated  and  that  the  total  reaction  energy  of  different  carbon  based  materials  can  be 
predicted  from  four  energy  parameters,  Esx,  Epx,  Egraphite,  and  local  relaxation  energy.  The 
parameters  Esx  and  Epx  depend  on  the  external  reactants  (x),  and  Egraphiie  and  local  relaxation 
energy  do  not  depend  on  the  reactants.  Therefore,  once  an  external  reactant  is  chosen,  we  can 
determine  Esx  and  Epx  from  DFT  calculations  of  its  reaction  with  the  graphite.  From  these  energy 
parameters,  we  can  predict  the  total  reaction  energies  and  binding  configurations  of  the  reactants 
on  diverse  fullemenes  and  nanotubes.  This  analysis  is  based  on  local  configurations  of  carbon 
atoms  so  that  one  can  apply  this  method  to  predict  local  deformation  effects  in  diverse  carbon 
nanostructres  without  performing  full  DFT  simulations. 
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Figure  3.  Hydrogenation  energies  of  graphene  sheet,  various  nanotubes,  and  C6o-  Circles  are 
DFT  simulation  data  with  full  relaxations,  and  squares  are  the  predicted  values  from  our  analysis. 
Data  points  represent  from  left  to  right,  graphene  sheet,  and  four  different  sizes  of  nanotubes 
listed  above,  and  Qo- 


ACKNOWLEDGMENTS 

This  work  is  supported  by  NASA  Ames  Director's  Discretionary  Fund  (DDF).  The  DFT 
simulations  are  performed  under  the  NPACI  allocation  SUA239  "Nanoscale  Materials 
Simulations." 


REFERENCES 

1 . H.  W.  Kroto,  J.  R.  Heath,  S.  C.  O’Brien,  R.  F.  Curl  and  R.  E.  Smally,  Nature  318,  162  (1985) 

2.  J.  Kong,  N.  R.  Franklin,  C.  W.  Zhou,  M.  G.  Chapline,  S.  Peng,  K.  Cho,  and  H.  Dai,  Science 
287,  622  (2000) 

3.  S.  J.  Tans,  A.  R.  M.  Verschueren,  and  C.  Dekker,  Nature  393,  49  (1998) 

4.  R.  C.  Haddan,  Science  261,  1545  (1993) 

5.  A.  Hirsch,  Topics  in  Current  Chemistry  198,  1 (1998) 

6.  M.  C.  Payne,  M.  P.  Teter,  D.  C.  Allan,  T.  A.  Arias  and  J.  D.  Joannopoulos,  Rev.  Mod.  Phys. 
64,1045(1992) 

7.  S.  Ismail-Beigi  and  T.  A.  Arias,  Comp.  Phys.  Comm.  128,  1 (2000) 

8.  R.  C.  Haddon,  Chem.  Phys.  Lett.  125,  231  (1986) 


Wl.5.6 


